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Abstract: Genetic variability and divergence among natural populations of Jaculus jaculus and J. orientalis in 
Tunisia were examined by electrophoretic analysis of 16 enzymatic proteins encoded by 23 genetic loci. Low levels of 
genetic variability were found among populations of both species in comparison to those of other rodent and mammal 
species of which data are available. In J. jaculus populations, the mean level of observed heterozygosity (Hops) ranged 
from 0.08 to 0.19, while the mean percentage of polymorphic loci (P) ranged from 26.2% to 45.2% and the mean number 
of alleles per locus (A) ranged from 1.1 to 1.4. Nevertheless, the mean values were 0.10 to 0.15, 29.3% to 44.1% and 1.1 
to 1.7, respectively, for J. orientalis. In addition, populations of the two species have revealed a lower degree of genetic 
differentiation (F,=0.0017 for J. jaculus and 0.0019 for J. orientalis). Moreover, F,, was 0.607, P<0.05 between 
populations of the two species, indicating that they were highly genetically diverged from each other. The present data 
assures the previous results on the validity of the present taxonomic situation of the two species and emphasis on the 
effect of geographic factors (environments type and bioclimatic stages) on the genetic structure of both species. 
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FRE: H 16 PE A Amt) 23 MEAR BE DT R Je AE Jaculus jaculus) MUR RHE. orientalis) AA 
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The genus Jaculus (Erxleben, 1777) basically 
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Israel and Morocco (Osborn & Helmy, 1980; Aulagnier 


contains two morphologically distinct species, lesser 
jerboa J. jaculus (Linnaeus, 1758) and greater jerboa J. 
orientalis (Erxleben, 1777), that occur in very diverse 
habitats throughout the sub-Saharan and deserts of North 
Africa, Asia, and Arabian countries such as Egypt, Sudan, 
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& Thévenot, 1986; Brown, 1994; Kingdon, 1997). The 
taxonomy of these species has received a considerable 
attention and has been the subject to controversial 
discussion. Pocock (1922) placed J. jaculus into the 
genus Jaculus and J. orientalis into the genus Scirtopoda 
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(Fischer, 
genitalia. However, on the basis of cranial and dental 


1817) due to differences in the external 


characters, Vinogradov (1930) reclassified J. orientalis to 
the genus Jaculus. This taxonomic controversy has 
evoked several additional studies on the Egyptian 
dipodids to clarify their taxonomic relationship. For 
example, osteology studies revealed that both J. jaculus 
and J. orientalis are more closely related to each other 
(Wassif, 1960). Chromosomal studies supported this 
hypothesis and interpreted the dissimilarity in the 
G-bands of the morphologically different chromosome 
pairs between the two species as pericentric inversions 
(Ata & Shahin, 1999; Ata et al, 2001; Shahin & Ata, 
2001). Moreover, 
deduced from allozymic survey (Shahin, 2003) have 


estimates of genetic divergence 


shown that J. orientalis appears to have shared a more 
recent common ancestor with J. jaculus than Allactaga 
tetradactyla. Divergence of these species would have 
occurred by Miocene (ca 9.6 to 18.7 million years ago). 
In Tunisia, the lesser jerboa J. jaculus is almost 
common in desert areas, while the greater jerboa J. 
orientalis occurs almost in the semiarid regions in the 
North. Due to the differences in the morphological 
features and ecological habitats of the two species as 
well as the lack of genetic investigations on both of them 
in Tunisia, the present study was undertaken to examine 


Zoological Research 


Vol. 30 


the allozymic variability and genetic divergence among 
populations of J. jaculus and J. orientalis. 


1 Materials and Methods 


1.1 Samples 

Adult individuals of both the lesser jerboa Jaculus 
jaculus and greater jerboa J. orientalis were collected in 
Tunisia between 2005 and 2007 from all currently known 
localities of their distribution as described by Burhan 
(1997). A total of 300 specimens collected from 12 
localities (six populations for each species) were live 
trapped and examined. Samples of J. jaculus were 
collected from Matmata, Nefta, Tataouine, Hamma, 
Gabes and Mednine, while those of J. orientalis were 
trapped from Mateur, El Khouat, Amra, Boumerdes, 
Oueslatia and Lessouda (Fig. 1). 
1.2 Methods 

Tissues from each specimen were preserved at 
—80°C until processed. Homogenates for electrophoresis 
were obtained from fractions of liver, muscle or kidney 
tissue crushed in distilled water. Horizontal starch gel 
(12%) electrophoresis of allozymes was carried out 
according to the protocols described in Harris & 
Hopkinson (1976) and Pasteur et al (1986). Twenty-three 
putative loci encoding 16 enzymatic proteins were 
analyzed using four different buffers (Tab. 1). 


Tab. 1 Enzymatic and non-enzymatic proteins surveyed among the populations of 
Jaculus jaculus and J. orientalis examined 











Enzyme Locus E.C Polymorphism Tissue Buffer, pH 

Isocitrate deshydrogenase ne 1.1.1.42 : Kidney TC, 6.7 

Malic enzyme Mod-1 1.1.1.40 P Kidney TC, 8.0 

Superoxide dismutase Sod-1 1.15.1.1 P Kidney TC, 6.7 

Glucose phosphate isomerase Gpi-1 5.3.1.9 P Kidney TC, 6..7 

Phosphogluconate dehydrogenase Pgd-1 1.1.1.44 P Liver TC, 8.0 

. Aat-1 P . 

Aspartate aminotransferase 2.6.1.1 Liver TME, 6.9 
Aat-2 P 

Glucose-6-phosphate dehydrogenase G6pd-1 1.1.1.49 P Liver TME, 6.9 

Sorbitol dehydrogenase Sdh-1 1.1.1.14 P Kidney TC, 8.0 

Alcohol dehydrogenase Adh-1 1.1.1.1 P Liver TME, 6.9 

Mannose phosphate isomerase Mpi-1 5.3.1.8 P Kidney TC, 6.7 

Esterase Es-2 3.1.1.1 P Kidney LI OH, 8.3 

Phosphoglucomutase Pgm-1 2.5.7.1 P Liver TME, 6.9 
Ldh-1 M 

Lactate dehydrogenase Ldh-2 1.1.1.27 M Kidney TC, 6.7 
Ldh-3 M 

Malate dehydrogenase ie late ip Kidney TC, 6.7 

Adenylate kinase ps RTAS n Kidney TC, 6.7 

Creatine kinase ar 2732 y Kidney TC, 6.7 
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Tab. 2 Allele frequencies at the 14 polymorphic loci surveyed among populations of Jaculus jaculus and 


J. orientalis examined 


Locality and sample size (N) 











J. orientalis J. jaculus 
El Bou- 

ae Mateur Khou- merd- ee Amra oe ori Nefta aa ee si Gabes 

NTA p Nes M25 NAS yeas; Neos: M35 yeas neas nens MP 
Est-2 
100 0.9750 0.7500 0.8000 0.7750 0.5750 0.5500 0.4000 0.4000 0.3750 0.3500 0.1750 0.3000 
110 0.0250 0.2500 0.2000 0.1750 0.4250 0.4000 0.1250 0.0250 0.0500 0.0250 0.0750 0.0000 
120 0.0000 0.0000 0.0000 0.0500 0.0000 0.0500 0.4750 0.5750 0.5750 0.6250 0.7500 0.7000 
Moa-1 
070 0.8250 0.8000 0.7000 0.7500 0.8250 0.5500 0.3250 0.3500 0.3250 0.3500 0.3000 0.1750 
080 0.1750 0.2000 0.2750 0.2500 0.1750 0.4500 0.0250 0.0000 0.1500 0.0250 0.0000 0.0750 
090 0.0000 0.0000 0.0250 0.0000 0.0000 0.0000 0.6500 0.6500 0.5250 0.6250 0.7000 0.7500 
Sod-1 
100 0.9250 0.7250 0.7500 0.8000 0.8500 0.9250 0.3500 0.3500 0.2500 0.2250 0.2000 0.0250 
110 0.0750 0.2750 0.2500 0.2000 0.1500 0.0750 0.6500 0.6500 0.7500 0.7750 0.8000 0.9750 
Idh-1 
100 0.9250 0.9750 0.8750 0.9500 0.8750 0.8500 0.2250 0.1500 0.2000 0.4250 0.1500 0.2250 
110 0.0750 0.0250 0.1250 0.0500 0.1250 0.1500 0.7750 0.8500 0.8000 0.5750 0.8500 0.7750 
Idh-2 
080 0.8250 0.9000 0.9250 0.8250 0.8750 0.9250 0.1750 0.1500 0.1250 0.1750 0.3000 0.2250 
090 0.1750 0.1000 0.0750 0.1750 0.1250 0.0750 0.8250 0.8500 0.8750 0.8250 0.7000 0.7750 
Pgd-1 
080 0.8250 0.5500 0.7250 0.7250 0.8000 0.7500 0.2750 0.1250 0.3000 0.1750 0.3500 0.3750 
090 0.1750 0.4500 0.2750 0.2750 0.2000 0.2500 0.0250 0.1000 0.0000 0.0750 0.0250 0.0500 
100 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.7000 0.7750 0.7000 0.7500 0.6250 0.5750 
Aat-1 
080 0.9000 0.8250 0.8250 0.8500 0.7750 0.8750 0.1750 0.1500 0.2250 0.1500 0.0500 0.8500 
090 0.1000 0.1750 0.1750 0.1500 0.2250 0.1250 0.8250 0.8500 0.7750 0.8500 0.9500 0.1500 
Aat-2 
100 0.8750 0.8500 0.8500 0.9250 0.8750 0.1250 0.2000 0.3000 0.0250 0.0500 0.0250 0.2250 
110 0.1250 0.1500 0.1500 0.0750 0.1250 0.8750 0.8000 0.7000 0.9750 0.9500 0.9750 0.7750 
G6pd1 
100 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
110 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 
Sdh-1 
090 0.4250 0.5750 0.5500 0.4750 0.8000 0.7750 0.1000 0.3000 0.0750 0.1750 0.1250 0.3500 
100 0.5750 0.3250 0.4500 0.4250 0.2000 0.1750 0.1000 0.0500 0.1250 0.1750 0.1000 0.0250 
110 0.0000 0.1000 0.0000 0.1000 0.0000 0.0500 0.4750 0.3250 0.5500 0.3250 0.7750 0.6250 
120 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.3250 0.3250 0.2500 0.3250 0.0000 0.0000 
Adh-1 
090 0.5750 0.5500 0.6500 0.5250 0.5500 0.4250 0.0500 0.1500 0.1250 0.0750 0.0750 0.1750 
100 0.4250 0.4000 0.2500 0.4000 0.4500 0.5750 0.0000 0.1000 0.0000 0.0500 0.1250 0.1750 
110 0.0000 0.0500 0.1000 0.0750 0.0000 0.0000 0.8000 0.5000 0.5000 0.5500 0.5500 0.3250 
120 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.1500 0.2500 0.3750 0.3250 0.2500 0.3250 
Mpi-1 
100 0.9250 0.8750 0.9000 0.7750 0.8250 0.7250 0.1250 0.0250 0.1500 0.1750 0.2250 0.0500 
110 0.0750 0.1250 0.1000 0.2250 0.1750 0.2750 0.8750 0.9750 0.8500 0.8250 0.7750 0.9500 
Pgm-1 
100 0.8750 0.8500 0.8500 0.9250 0.8750 0.2250 0.2000 0.3000 0.0250 0.1250 0.0250 0.0500 
110 0.1250 0.1500 0.1500 0.0750 0.1250 0.7750 0.8000 0.7000 0.9750 0.8750 0.9750 0.9500 
Gpi-1 
110 0.9250 0.9000 0.1250 0.9250 0.8750 0.8500 0.8500 0.2250 0.8000 0.6500 0.2250 0.8500 
120 0.0750 0.1000 0.8750 0.0750 0.1250 0.1500 0.1500 0.7750 0.2000 0.3500 0.7750 0.1500 


1.3 Data analysis 


Observed electromorphs or alleles were identified 
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herein according to their electrophoretic mobility relative 


to that of the most common electromorph (assigned 
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mobility=100). Gene frequencies and genetic variability 
parameters were calculated using the GENETIX 4.03 
software (Belkhir et al, 2001). The standardized genetic 
variance (F«) among all populations was estimated for 
the polymorphic loci and the significance of the Fs 
values was tested using Monte Carlo simulations in 
ARLEQUIN 3.0 (Excoffier et al, 2005). Nei’s (1972) 
genetic distance (D) matrices between populations were 
8°E 


$ Mateu 
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calculated using the program GENDIST from the 
PHYLIP 3.5 package (Felsenstein, 1993), which was 
then used for the construction of the phenogram by the 
UPGMA (Sneath & Sokal, 1973). Bootstrap values were 
obtained from 1000 pseudo replicates of allele 
frequencies using the SEQBOOT routine in PHYLIP. 
The sequential Bonferroni’s test (Rice, 1989) was used to 
correct for multiple tests. 
10°E 12°E 
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Fig. 1 A map showing the geographical localities from which the populations of both Jaculus jaculus 


and J. orientalis were collected in Tunisia 


2 Results 


Twenty-three loci encoding 16 enzymatic systems 
were compared among populations of the two studied 


jerboa species. Of these 23 loci, only nine (39%) loci 
(Ldh-1, Ldh-2, Ldh-3, Mdh-1, Mdh-2, Ak-1, Ak-2, Ck-1 
and Ck-2) were monomorphic with the same allele fixed 
in all populations, while the remaining 14 (61%) loci 
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were polymorphic with different alleles (Tab. 2). 
2.1 Intraspecific genetic structure 

A summary of genetic data on the 12 populations of 
the two dipodid species is given in Tab. 3. The observed 
heterozygosity (H.»;) ranged from 0.10 to 0.15 in J. 
jaculus, while it varied between 0.08 and 0.19 in J. 
orientalis. However, larger values of the percentage of 
polymorphic loci (P) were noted in both species; it 
ranged from 27.28 % to 44.11 % in J. orientalis and 
22.12 % to 45.19 % in J. jaculus. Similarly, allelic 
richness (number of alleles per locus (A) exhibited 
higher levels in both J. jaculus (range from 1.09 to 1.65) 
and J. orientalis (range from 1.11 to 1.42). 

It is worthy to mention that an absence of the 
genetic differentiation was noted in the populations of 
the two species. This was clearly demonstrated by the 
relatively lower levels of Fy, and Nei’s distance. In the 
lesser Jerboa J. jaculus, F value was 0.0017 and Nei’s 
distance (0.0016), while in the greater jerboa J. orientalis, 
they were 0.0019 and 0.0018, respectively. These lower 
values of both Fy and Nei’s distance in both species 
reflect a high genetic homogeneity and a clear high gene 
flow among their populations (Tab. 4). 

2.2 Interspecific genetic relationship 

A clear difference in allelic frequencies and 
distribution was noted between the two species. In 
addition, the two species were genetically highly 
differentiated and the Fy, and Neis (1972) distance 
between the two species were 0.61 and 0.56, respectively. 
A diagnostic locus (G6pd-1) as well as remarkable 
differences in allelic frequencies and distribution were 
noted between the two species and yielded two separate 


clades by means of UPGMA analysis corresponding to 
the lesser jerboa and greater jerboa, respectively (Fig. 2). 


3 Discussion 


Allozyme analysis of the 23 genetic loci showed 
that the values of genetic variation (Hos, P and A) 
observed in Jaculus jaculus and J. orientalis are 
generally within the range reported for other rodent 
species of which data are available (Nevo et al, 1990) 
and slightly higher than those reported for the same 
species in Egypt (Shahin, 2003). 

Analysis of our data suggests that only low 
differentiation were found between local populations of 
the two species, as exemplified by the low values of 
genetic distances and F«. This result of close similarity in 
genetic content between local populations of the same 
taxon seemingly characterizes all types of organisms for 
which further data are available (Selander & Johnson, 
1973; Avise, 1974; Ayala, 1975). Estimates of levels of 
genetic similarity suggest that the dipodids examined 
here were comparatively either within the range or quite 
different from other taxa (Nevo et al, 1974; Avise & 
Smith, 1977; Gardenal et al, 1990; De Sousa et al, 1996). 
Significant 


correlation between levels of genetic 


diversity and ecological parameters (life zone, 
geographical range, habitat type and range, and climate 
region), demographic parameters (species size and 
population structure, gene flow and sociability), and a 
series of life history characteristics (longevity, generation 
length, fecundity, origin and parameters related to the 
mating system mode of reproduction) has previously 


been reported (Nevo et al, 1984). 


Tab.3 A summary of genetic variability in the populations of Jaculus jaculus and J. orientalis 
examined based on electrophoretic analysis of 23 loci 




















Mean no. of Percentage of Mean heterozygosity 
Species Populations alleles per polymorphic loci 
locus (A) (P) Hovs Hop 
Mateur 1.23 44.11 0.1114 + 0.0378 0.1371 + 0.0280 
El Khouat 1.17 34.15 0.0950 + 0.0112 0.1244 + 0.0528 
Boumerdes 1.43 29.27 0.1171 + 0.0646 0.1140 + 0.0266 
J. orientalis Oueslatia 1.24 41.18 0.0964 + 0.0646 0.1173 + 0.0426 
Amra 1.65 35.98 0.1536 + 0.0308 0.1884 + 0.0237 
Lessouda 1.09 27.28 0.1314 + 0.0469 0.1707 + 0.0695 
Mean 1.31 35.33 0.1175 + 0.0427 0.1420 + 0.0405 
Matmata 1.42 31.18 0.1629 + 0.0132 0.1996 + 0.0730 
Nefta 1.11 45.19 0.1886 + 0.0145 0.1927 + 0.0550 
Tataouine 1.23 22.12 0.0814 + 0.0469 0.1084 + 0.0237 
J. Jaculus Mednine 1.39 28.15 0.1929 + 0.0132 0.2007 + 0.0695 
Hamma 1.16 26.18 0.1386 + 0.0145 0.1896 + 0.0730 
Gabes 1.29 44.01 0.1186 + 0.0145 0.1727 + 0.0550 
Mean 1.27 32.81 0.1472 + 0.0195 0.1773 + 0.0582 
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Tab. 4 Pairwise estimates of F „values between the populations of Jaculus jaculus and J. orientalis based on 
the 14 polymorphic loci are given above (diagonal) and the genetic distances (Nei, 1972) are below 
(diagonal). The correct values of sequential Bonferroni’s test between the two species are added 


J. orientalis 




















Populations 
Mateur El Khouat Boumerdes Oueslatia Amra Lessouda 
Mateur = 0.0014 0.0028 0.0016 0.001 0.0013 
El Khouat 0.0025 S 0.0024 0.0029 0.0012 0.0017 
T orientalis Poum 0.0027 0.0014 = 0.0012 0.0025 0.0026 
Oueslatia 0.0019 0.0012 0.0021 = 0.0027 0.0023 
Amra 0.0028 0.0019 0.0014 0.0023 = 0.0011 
Lessouda 0.001 0.0022 0.0011 0.0016 0.0023 = 
Matmata 0.5125* 0.6116* 0.4567* 0.5431* 0.5987* 0.6532* 
Nefta 0.5234* 0.5678* 0.6433* 0.4765* 0.5378* 0.6224* 
re Tataouine 0.5532* 0.4987* 0.5342* 0.5643* 0.5899* 0.4568* 
Mednine 0.4912* 0.5112* 0.6312* 0.6435* 0.5668* 0.6345* 
Hamma 0.4669* 0.5678* 0.5631* 0.6987* 0.4889* 0.5766* 
Gabes 0.5678* 0.4876* 0.6111* 0.5467* 0.5346* 0.5128* 
J. jaculus 
Populations 
Matmata Nefta Tataouine Mednine Hamma Gabes 
Mateur 0.3356* 0.4567* 0.5676* 0.6570* 0.7840* 0.6573* 
El Khouat 0.6634* 0.5637* 0.6543* 0.5768* 0.7654* 0.5678* 
Gt eae Boumerdes 0.5564* 0.4967* 0.3456* 0.7890* 0.5678* 0.6578* 
Oueslatia 0.6546* 0.6653* 0.6748* 0.5678* 0.4567* 0.5578* 
Amra 0.4563* 0.7654* 0.5678* 0.6543* 0.6784* 0.6789* 
Lessouda 0.4650* 0.8754* 0.4567* 0.6789* 0.5671* 0.7689* 
Matmata -= 0.0012 0.0016 0.0019 0.0027 0.0022 
Nefta 0.0012 -= 0.0001 0.0028 0.0008 0.0021 
. Tataouine 0.0018 0.0014 = 0.0024 0.0011 0.0024 
J. jaculus : 
Mednine 0.001 0.002 0.0013 = 0.0013 0.0017 
Hamma 0.0021 0.0014 0.0024 0.0027 = 0.0029 
Gabes 0.0013 0.0011 0.0026 0.0022 0.0031 0.0016 


* Significant after correction Bonferronni’s. 


In the present investigation, similar explanation 
linked to genetic variability and genetic differentiation in 
J jaculus and J. orientalis, could be assumed when 
comparing the present data with that of the Egyptian 
species (Shahin, 2003). 

The main causes of genetic variability and relatively 
low heterozygosity exhibited by the studied populations 
have been discussed by several authors. For example, it 
has been reported that both random and deterministic 
factors, including genetic drift, selection, migration, 
mutation and historic events, may affect the population 
size and breeding and thereby causing homozygosity and 
reducing heterozygosity (Nevo et al, 1974). In addition, 
an increase in genetic variability could be adaptative 
strategy in an unexpected environment (Nevo, 1978), and 
variability could also remain weak in an ecologically 
diversified environment (Pasteur et al, 1978). Moreover, 
under stable conditions in a uniform trophic environment, 
genetic variability could accumulate (Ayala & Valentine, 
1974). Furthermore, it has been pointed out that the 
genetic variation observed among populations living in 


nearly stable environmental conditions could be 
suggestive of differences in vagility and breeding 
(Gorman et al, 1977). Normally, high vagility and 
consequent low inbreeding results in relatively high 
levels of genetic variation. Thus, like in many other 
rodent species, the considerable levels of genetic 
variability observed in these dipodid species could be 
explained as an adaptative strategy for homozygosity in 
the relatively uniform environment (Shahin, 2003). On 
the other hand, the significant increase of heterozygotes 
more than expected among the dipodid species examined 
in this study may be due to: selection, either for the 
heterozygotes or against the homozygotes, negative 
assortative mating, or any other special explanation, that 
the heterozygotes are more active than homozygotes and 
that they are frequently more trapped (Shahin, 2003). 
Nevertheless, the low intraspecific Fs and Nei’s distance, 
observed in this study, may be explained by high gene 
flow among populations linked probably to high vagility 
of both species (Anderson, 1970; Wolff, 2008). These 
findings are in accordance with those found in Egyptian 
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Fig.2 Aphenogram constructed based on Nei's (1972) genetic distance (D) matrices between 
the populations of Jaculus jaculus and J. orientalis examined 


counterparts (Shahin, 2003). When 


distances are compared between species, higher values 


interspecific 


were noted indicating probably the effect of the 
bioclimatic factors on the genetic differentiation between 
species. Indeed, J. jaculus occupies the southern parts of 
Tunisia, mainly desert areas, while J. orientalis is 
common in the north of Tunisia occupying mountains 
and dense vegetations. 

In conclusion, the genetic relationship between the 
examined, as 


two species demonstrated by the 


phylogenetic tree, is generally similar to that reported for 


the same species in Egypt (Shahin, 2003) and indicates 
that the divergence of the two species from their 
common ancestor has occurred since 3 millions years 
ago. In addition, the present data assures the validity of 
the present taxonomic situation of the two species and 
effect 
(environments type and bioclimatic stages) on the 


emphasis on the of geographic factors 


genetic structure of both species. 
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